Non-adiabatic Effects in the Braiding of Non-Abelian Anyons in Topological 

Superconductors 
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Qubits in topological quantum computation are built from non-Abelian anyons. Adiabatic braid- 
ing of anyons is exploited as topologically protected logical gate operations. Thus, the adiabatic- 
ity upon which the notion of quantum statistics is defined, plays a fundamental role in defining 
the non-Abelian anyons. We study the non-adiabatic effects in braidings of fsing-type anyons, 
namely Majorana fermions in topological superconductors, using the formalism of time-dependent 
Bogoliubov-de Gennes equations. Using this formalism, we consider non-adiabatic corrections to 
non-Abelian statistics from: (1) tunneling splitting of anyons imposing an additional dynamical 
phase to the transformation of ground states; (2) transitions to excited states that are potentially 
destructive to non-Abelian statistics since the non-local fermion occupation can be spoiled by such 
processes. However, if the bound states are localized and being braided together with the anyons, 
non-Abelian statistics can be recovered once the definition of Majorana operators is appropriately 
generalized taking into account the fermion parity in these states. On the other hand, if the excited 
states are extended over the whole system and form a continuum, the notion of local fermion parity 
no longer holds. We then quantitatively characterize the errors introduced in this situation. 

PACS numbers: 03.65.Lx, 05.30.Pr, 03.67.Lx, 03.67.Pp 



I. INTRODUCTION 



Fragility of quantum coherence against external per- 
turbations is one of the main obstacles in building a quan- 
tum computer. Maintaining coherence in the presence 
of environment is thus among the very first challenges 
in the realization of quantum computation. Topological 
Quantum Computation(TQC)i2 provides an elegant so- 
lution to the problem: quantum information is encoded 
and manipulated in a non-local way such that weak local 
perturbation has essentially no influence on the compu- 
tation procedure so the decoherence is eliminated at the 
hardware level. To achieve this goal, qubits are built out 
of non-Abelian anyons, which are exotic particles obey- 
ing non-Abelian statistics, and adiabatic braidings of the 
anyons generate quantum gates. The intrinsic protec- 
tion of quantum coherence against environmental noise 
makes TQC a very appealing scheme for quantum com- 
putation. It is thus very desirable to understand the 
properties of non-Abelian anyons as well as to search for 
quantum phases of matter that host them, known as non- 
Abelian topological phases. 

Quantum statistics describes how many-body wave- 
functions transform under the exchange of any pair of 
particles. The unitary transformations form represen- 
tations of the braid group^ In (2+l)-dimensions, non- 
Abelian statistics happens when there are degenerate 
ground states, and braidings (pair- wise exchanges) result 
in unitary rotations in the ground state subspace, cor- 
responding to higher-dimensional representation of the 
braid group in two-dimensions. These braiding opera- 
tions form the basic set of topologically protected gates 



to perform quantum computation. The simplest of such 
anyons, and probably the only one which can be po- 
tentially realized, is the so-called Ising-type anyons. It 
has a natural interpretation as Majorana fermions which 
are self-conjugate neutral fermionic excitations. Such 
excitations are found to exist in some topological su- 
perconducting systems such as two-dimensional spin- 
less chiral p-wave superconductors. They are in one- 
to-one correspondence with non-degenerate zero-energy 
states in superconductors which exist in vortices in two- 
dimensionsi— or domain walls in one-dimension?. Non- 
Abelian statistics of Majorana fermions in both cases 
have been demonstrated explicitly 7 -^ and are also sup- 
ported by rather general arguments based on the con- 
servation of fermion parityii. For experimental real- 
izations, it has long been conjectured that quasiholes 
in v = 5/2 fractional quantum Hall state are Ising 
anyon a 12 " 14 and subsequently candidate materials such 
as Sr2Ru04 for chiral p-wave superconductivity were also 
discovered^. However, recent theoretical advances have 
provided a variety of different heterostructure systems, 
where spin-orbit coupled semiconductors^— or the sur- 
faces of three-dimensional topological insulators^ com- 
bined with proximity-induced s-wave pairing are shown 
to successfully engineer spinless p x + ip y superconductiv- 
ity with non-Abelian Majorana anyons. 

Mathematical definition of quantum statistics neces- 
sarily builds upon the concept of Berry phase of many- 
body wavefunctions. This implies that the adiabaticity of 
braiding is an essential ingredient for non-Abelian statis- 
tics, since the quantum state has to stay in the ground 
state manifold during the entire process of the braiding. 
In the real world, however, braidings are necessarily per- 
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formed within a finite time interval, i.e., they are al- 
ways non-Adiabatic. As known from the adiabatic per- 
turbation theory, Berry phase is the leading-order term 
in the adiabatic perturbative expansion^ - — Given the 
fundamental role played by adiabatic braiding in TQC, it 
is therefore important to understand quantitatively the 
higher order corrections arising from non-adiabatic evo- 
lution, 

In this paper, we present a systematic study of the 
non-adiabatic corrections to the braiding of non-Abelian 
anyons and develop formalism to describe their dynam- 
ical aspects. In our treatment, braidings are considered 
as dynamical evolutions of the many-body system, es- 
sentially using the time-dependent Schrodinger equation 
of the BCS condensate whose solutions are derived from 
time-dependent Bogoliubov-de Gennes (BdG) equation. 
Generally, adiabaticity may break down in three differ- 
ent ways: (a) tunneling of non-Abelian anyons when 
there are multiples of them, which splits the degener- 
ate ground state manifold and therefore introduces ad- 
ditional dynamical phases in the evolution; (b) transi- 
tions to excited bound states outside the Hilbert space 
of zero-energy states, in this case topologically protected 
braidings have to be defined within an enlarged Hilbert 
space; (c) transitions to the continuum of extended states 
which render the fermion parity in the low-energy Hilbert 
space ill-defined. These non-adiabatic effects are possible 
sources of errors for quantum gates in TQC. The main 
goal of this paper is to quantitatively address these effects 
and their implications on quantum computation. 

Our work is the first systematic attempt to study 
non-adiabaticity in the anyonic braiding of non-Abelian 
quantum systems. Given that the braiding of non- 
Abelian anyons is the unitary gate operation— in topo- 
logical quantum computation--, understanding the dy- 
namics of braiding as developed in this work is one of 
the keys to understanding possible errors in topologi- 
cal quantum computation. The other possible source 
of error in topological quantum computation is the lift- 
ing of the ground state anyonic degeneracy due to inter- 
anyon tunneling, which we have studied elsewher e 25 ' 26 . 
Although we study the braiding non-adiabaticity in the 
specific context of the topological chiral p-wave super- 
conductors using the dynamical BdG equatons within 
the BCS theory, our work should be of general va- 
lidity to all known topological quantum computation 
platforms, since all currently known non-Abelian any- 
onic platforms in nature are based on the SU(2)2 con- 
formal field theory of Ising anyons, which are all iso- 
morphic to the chiral p-wave topological superconduc- 
tors 2 -. As such, our work, with perhaps some mi- 
nor modifications in the details, should apply to the 
fractional quantum Hall non-Abelian qubits 2 ^, real p- 
wave superconducting systems based on solids 2 - and 
quantum gasea 2 ^, topological insulator-superconductor 
heterostructures 2 ^, and semiconductor-superconductor 
sandwich structures^ 6 - and nanowire a 18 ' 29 . Our results 
are quite general and are independent, in principle, of the 



detailed methods for the anyonic braiding which could 
vary from system to system in details. 

The reminder of the paper is organized as follows: in 
Sec. [TT] we review Majorana fermions in non-Abelian 
topological superconductors and their quantum statis- 
tics defined in terms of adiabatic evolution. In Sec. IIIII 
we generalize BdG equation to describe time-evolution 
of BCS superconductors under a parametrically time- 
dependent Hamiltonian. Then we study the three main 
non-adiabatic effects as outlined above using the ap- 
proach of time-dependent BdG equation. 



II. QUANTUM STATISTICS OF MAJORANA 
FERMIONS 

A. Quantum Statistics and Adiabatic Evolution 

We first briefly review how quantum statistics is for- 
mulated mathematically in terms of the adiabatic evo- 
lution of many-body wavef unctions, following a recent 
exposition in Ref. [14j . Consider the general many-body 
Hamiltonian W [Ri (i) , . . . , R„(i)] where parameters {R;} 
represent positions of quasiparticles. We assume the ex- 
istence of well-defined, localized excitations which we call 
quasiparticles. At each moment t, there exists a subspace 
of instantaneous eigenstates of H[R.i(t), . . . , R„(i)] with 
degenerate energy eigenvalues. Instantaneous eigenstates 
in the subspace are labeled as \a(t)) = |a({R; (£)})). We 
constraint our discussion in the ground state subspace 
with zero energy eigenvalue. 

The adiabatic exchange of any two particles can be 
mathematically implemented by adiabatically changing 
the positions of two particles, say R; and Rj, in such a 
way that in the end they are interchanged. This means 
that 

R i (T) = R i (0),R,(T) = R i (0). (1) 

According to the adiabatic theorem^ ., it results in a uni- 
tary transformation within the subspace: if the system 
is initially in state |^(0)), then \tp{T)) = U\i/)(0)). To 
determine U, we first consider initial states |-0(O)) = 
|a({R,(0)})}. Under this evolution the final state can 
be written as 

\MT)) = UoWT)). (2) 

Here the matrix Uq is the non-Abelian Berry 
phase; 21 ' 31 ^ 2 - 

U = Venp^ij^ dtM(t)j , (3) 

where V denotes path ordering and matrix element of the 
Berry's connection M is given by 

M afj (t)=i{a(t)\p(t)). (4) 
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Although the exchange defines a cyclic trajectory in 
the parameter space of Hamiltonian, the final basis states 
can be different from the initial ones (e.g, |a({Rj})) can 
be multivalued functions of R^, which is allowed if we are 
considering quasiparticles being collective excitations of 
many-body systems). The only requirement we impose 
is that the instantaneous eigenstates {a(t)} are contin- 
uous in t. Therefore, we have another matrix B de- 
fined as B afi = (a(0)\P(T)), relating {a(T)} to {a(0)}: 
\a(T)) = B a p\j3(0)). Combining with ([2]), we now have 
the expression for U: 

\i>(T)) = U B^(0)). (5) 

Therefore 

U = U B = Pexp^i J dtM(t)^j B. (6) 

In fact, the factorization of U into Uq and B is somewhat 
arbitrary and gauge- dependent. However, their combi- 
nation U is gauge-independent provided that the time- 
evolution is cyclic in parameter space (positions of parti- 
cles). The unitary transformation U defines the statistics 
of quasiparticles. 

B. Non-Abelian Majorana Fermions 

We now specialize to the non-Abelian statistics of Ma- 
jorana fermions in topological superconductors, carefully 
treating the effect of Berry phases. We mainly use spin- 
less superconducting fermions as examples of topological 
superconductors in both ID and 2D since all known topo- 
logical superconducting systems supporting non-Abelian 
excitations essentially stem from spinless chiral p-wave 
superconductors i 6 ' 17 i 20 

In the BCS mean-field description of superconduc- 
tors, the Hamiltonian is particle-hole symmetric due to 
U(l) symmetry breaking. In terms of Nambu spinor 
i&(r) = (?/>(r), ip (r)) T , the BCS Hamiltonian is expressed 
as "Hbcs = I Jd 2 r^ t (r)iJ B dG*(r)- The Bogoliubov-de 
Gcnncs Hamiltonian -HBdG takes the following form£i22 

H BdG = ^ t _t T ) ' ( 7 ) 

where h is the single-particle Hamiltonian [for spinless 
fermions it is simply h — (-^9^ — /it) 5(r — r')] and 
A is the gap operator. The BCS Hamiltonian can be 
diagonalized by Bogoliubov transformation 

7 f = J d 2 r [u(r)v3 f (r) + u(r)^(r)] . (8) 

Here the wavefunctions u(r) and v(r) satisfy BdG equa- 
tions: 

M$D -*($)• (9 ) 



Throughout this work, we adopt the convention that op- 
erators which are hatted are those acting on many-body 
Fock states while bold ones denote matrices in "lattice" 
space. 

The single-particle excitations 7, known as Bogoliubov 
quasiparticles, are coherent superpositions of particles 
and holes. The particle-hole symmetry implies that the 
quasiparticle with eigenenergy E and that with eigenen- 
ergy —E are related by J-e = 1 E - Therefore, E = 
state corresponds to a Majorana fermion 70 = 7q 8 . The 
existence of such zero-energy excitations also implies a 
non-trivial degeneracy of ground states: when there are 
2A such Majorana fermions, they combine pair-wisely 
into N Dirac fermionic modes which can either be occu- 
pied or unoccupied, leading to 2 N -io\d degenerate ground 
states. The degeneracy is further reduced to 2 N ~ 1 by 
fermion parity^. Since these fermionic modes are intrin- 
sically non-local, any local perturbation can not affect 
the non-local occupancy and thus the ground state de- 
generacy is topologically protected. This non-locality lies 
at the heart of the idea of topological qubits. 

We are mostly interested in Majorana zero-energy 
states that are bound states at certain point defects {e.g. 
vortices in 2D, domain walls in ID). In fact, Majorana 
bound states are naturally hosted by defects because that 
zero-energy states only appear when gap vanishes. De- 
fects can be moved along with the Majorana fermions 
bound to them. Braidings of such Majorana fermions 
realizes very non-trivial non-Abelian statistics. 

We now apply the general theory of quantum statis- 
tics as previously discussed in Sec. Ill Al to the case of 
Majorana fermions in topological superconductors. The 
simplest setting where non-trivial statistics can be seen 
is the adiabatic braiding of two spatially separated Majo- 
rana fermions 71 and 72. We denote the two bound state 
solutions of the BdG equation by ^01 ancl *02- When 
Ri and R2 vary with time they become instantaneous 
zero-energy eigenstates of BdG Hamiltonian. We choose 
their phases in such a way that the explicit analytical 
continuation of BdG wavefunction leads to the following 
basis transformation under exchange 7 

*0l(T) = 5*02(0) 

* 02 (T) = -s* i(0), 1 ' 

where s = ±1. The value of s depends on the choice of 
wavefunctions and we choose the convention that s — 1 
throughout this work. In the case of Majorana fermions 
in vortices, the additional minus sign originates from 
branch cuts introduced to define the phase of wavefunc- 
tions. This transformation actually gives the B matrix in 
the general theory. Equivalently in terms of quasiparticle 
operator, we have 

* 72 
72 — > —71 ■ 

If we define the non-local fermionic mode d ' = (71 + 
172), the states with even and odd fermion parity are 
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given by |g) and cft\g). Due to the conservation of fermion 
parity, the two states are never coupled. However, the 
non-Abelian statistics still manifest itself in the phase 
factor acquired by the two states after an adiabatic ex- 
change. To see this, first we notice that under exchange, 
the analytical continuation(or basis transformation) gives 
the following transformation of the two states: 



->e»*|g) 



-> e 2 e 



(12) 



Here the ^ phase difference is reminiscence of non- 
Abelian statistics. 

So far we have obtained the basis transformation ma- 
trix B. To know the full quantum statistics we also need 
to calculate the adiabatic evolution Uq. We now show by 
explicit calculation that Uq cx 1 up to exponentially small 
corrections. This requires knowledge of Berry connec- 
tion accompanying adiabatic evolutions of BCS states. 
Fortunately, for BCS superconductor the calculation of 
many-body Berry phase can be done analytically 3 -^. The 
ground state |g) has the defining property that it is anni- 
hilated by all quasiparticle operators j n . All other states 
can be obtained by populating Bogoliubov quasiparticlcs 
on the ground state |g). Let us consider a state with M 



quasiparticles \ni, n 2 , . . . , n M ) = T ni T n2 • ■ • % M |g)- The 
Berry connection of this state then read: 



:S4 



(ni, ... , UM\d\nx, . . . , n M ) 



m , v 

(g\d\g)+Y,{<o<) d { v ni )- ( 13 ) 

So the difference between the Berry phase of a state with 
quasiparticles and the ground state is simply the sum of 
"Berry phase" of the corresponding BdG wavefunctions. 
Since the Berry phase of ground state |g) can be elimi- 
nated by a global U(l) transformation, only the difference 
has physical meaning. 

According to (p~3|) , the relevant term to be evaluated is 

"01 + * w 02 
V 1 + 1^02, 

2Re (u*dui + u*du 2 ) + 2iRe(uldu 2 - u* 2 dui), (14) 

where we have made use of the Majorana condition 
v = u* . The first term in (fl"4"|) vanishes because J u*du 
must be purely imaginary. The second term has a non- 
vanishing contribution to the total Berry phase. How- 
ever, due to the localized nature of zero-energy state, the 
overlap between u\ and u 2 is exponentially small: 



{u* QX -iu*Q 2 ,vl x -iv* m )d 



r 

I dtKc(uldtu 2 — u 2 dtU\) 
Jo 



= -|Ri-R 2 |/£ 



(15) 



Therefore the Berry phase can be neglected in the limit of 
large separation R. This completes our discussion of non- 
Abelian statistics. The above calculation can be easily 
generalized to the case of many anyons. 



III. TIME-DEPENDENT BOGOLIUBOV-DE 
GENNES EQUATION 

In this section we derive the formalism to track down 
time-evolution of BCS condensate within mean-field the- 
ory. For BCS superconductivity arising from interac- 
tions, the pairing order parameter has to be determined 
sclf-consistently, which makes the mathematical problem 
highly nonlinear. In the situations that we are interested 
in, it is not critical where the pairing comes from. In some 
systems that are believed to be experimentally accessible, 
e.g. semiconductor/superconductor heterostructure, su- 
perconductivity is induced by proximity effec t 16 1 20 and 
there is no need to keep track of the self-consistency. We 
will take the perspective that order parameter is simply 
a external field in the Hamiltonian. 

The time-dependent BdG equatio n 35 i 36 has been 
widely used to describe dynamical phenomena in BCS 
superconductors. To be self-contained here we present 
a derivation of the time-dependent BdG equation high- 
lighting its connection to quasiparticle operators. It can 
also be derived by methods of Heisenberg equation of 
motion or Green's function. Suppose we have a time- 
dependent BdG Hamiltonian HBdcit)- The unitary time- 
evolution of the many-body system is formally given by 



U(t) = T exp 



-t I di'HBCs(t') 
/o 



(16) 



To obtain an explicit form of U (t) , we define the time- 
dependent Bogoliubov quasiparticle operator as 



7„(t) = U{t)%il\t), 



(17) 



where 7„ is the quasiparticle operator for Hbcs(0) and 
the corresponding BdG wavefunction is u n (r) , v n (r) . We 
adopt the normalization condition 



d 2 rK(r)| 2 + K(r)| 2 



(18) 



which means {j n ,7n} = 1. 

■J n (t) by definition satisfies the following equation of 
motion 



1 ^ = [ftBCs0),7n(*)]< 



(19) 



In fact, by direct calculation one can show that the 
equation of motion is solved by 



dr 



u n (r, t)i>(r) + iv, (r, t)tp(r) 



(20) 



where the wavefunction u n (r, t) and v n (r, t) are solutions 
of time-dependent BdG equation: 



d / u„(r, t)\ , u n (r, i) N 



(21) 
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together with initial condition (another way of saying 
7n(0) = 7n) 

u n (r, 0) = un(r), w B (r, 0) = v n (r). (22) 

As long as the solutions of the time-dependent BdG equa- 
tion (|2ip are obtained, we can construct the operators 

mm. 

We now derive an explicit formula of U when the time- 
evolution is cyclic(ie., Hbcs{T) = %bcs(0))- In that 
case, it is always possible to express 7n(T) as a linear 
combination of j n = j n (0) and 7^ . Since particle number 
is not conserved, it is more convenient to work with Ma- 
jorana operators. We thus write j n = C2 n —x +*C2n where 
c m are Majorana operators. Suppose the BdG matrix 
has totally 2N eigenvectors so n = 1,2 ... ,N. Assume 
that by solving time-dependent BdG equation we obtain 
the transformation of c m as follows: 

Ck(T) = J2 V kiCh (23) 
1 

where V G SO(2A) as required by unitarity and the 
conservation of fermion parity. The matrix V can be 
calculated once we know the BdG wavef unctions. Then 
we can write down an explicit expression of U (T) in terms 

of CrJ-: 

U{T) = CXp I 1 D rnnC m C n I , (24) 
\ ran J 

where the matrix D is defined by the relation e~ D = V. 
Here D is necessarily a real, skew-symmetric matrix. A 
proof of this result is given in Appendix [A] We notice 
that the usefulness of is actually not limited to cyclic 
evolution. In fact, (|24| is purely an algebraic identity 
that shows any SO(2iV) rotation of 2N Majorana opera- 
tors can be implemented by a unitary transformation. 

In the following we outline the method to solve the 
time-dependent BdG equation. To make connection with 
the previous discussion of Berry phase, we work in the 
"instantaneous" eigenbasis of time-dependent Hamilto- 
nian. At each moment t, the BdG Hamiltonian if BdG (i) 
can be diagonalized yielding a set of orthonormal eigen- 
functions {^E'„(r, t)}. A remark is right in order: because 
of particle-hole symmetry, the spectrum of BdG Hamil- 
tonian is symmetric with respect to zero energy and the 
quasiparticle corresponding to negative energy are really 
"holes" of positive energy states. However, at the level 
of solving BdG equation mathematically, both positive 
and negative energy eigenstates have to be retained to 
form a complete basis. The most general form of BdG 
wavefunction can be expanded as 

*(r,t) = 5^c n (t)*„(r,t). (25) 

a 

Plugging into time-dependent BdG equation, we ob- 
tain 

iCn + ^2 M nm (t)c m = E n (t)c n , (26) 

m 



where M nm = i J d 2 r*t (r)* m (r). 

Assume that starting from initial condition c„(0) = 
$mn (roughly the quasiparticle is in the ^ m state at 
t — 0), we obtain the solutions of (|2l)|) at t = T denoted 
by c™(T). The transformation of basis states themselves 
is given by the matrix B. Combining these two transfor- 
mations we find 

J2 c kBki% (27) 
kl 

from which the linear transformation V can be directly 
read off. Then by taking the matrix log of V we can ob- 
tain the evolution operator. This is the procedure that 
we will use to solve the (cyclic) dynamics of BCS super- 
conductors. 

We will not be attempting to obtain the most gen- 
eral solution, since it depends heavily on the microscopic 
details. Instead, we focus on two major aspects of non- 
adiabaticity: (a) finite splitting of ground state degener- 
acy which only becomes appreciable when the braiding 
time is comparable to the "tunneling" time of Majorana 
fermions. (b) excited states outside the ground state sub- 
space. 

Since these effects are highly non-universal, we have 
to specify how the braiding is performed. For example, 
when considering the braiding of non-Abelian vortices 
in two-dimensions, we assume that the trajectory of the 
two vortices is a circle and the braiding is performed at 
a constant speed. Mathematically, we choose Ri(t) and 
R-2(t) as 

Ri(t) = -R 2 (t) = R{cos(ujt + 6 ), sinM + 6 )). (28) 

Here we have defined to = I? characterizing the energy 
scale corresponding to non-adiabaticity of the braiding 
itself. 



A. Effect of tunneling splitting 

The derivation of transformation rule (jlip assumes 
that the two Majorana bound states have vanishing 
energies so there is no dynamical phase accumulated. 
The assumption is only true when tunneling splitting 
of zero-energy states is neglected. It has been estab- 
lished that the finite separation between anyons always 
leads to a non-zero splitting of zero-energy states£ ' 25 i 37 
although the splitting is exponentially suppressed due 
to the existence of bulk gap. As a result, the two 
ground states acquire different dynamical phases dur- 
ing the time-evolution. Here we take into account all 
the non-universal microscopic physics including dynami- 
cal phase induced by tunneling splitting and non-Abelian 
Berry phase. 

In the framework of time-dependent BdG equation, the 
two basis states are 1 $ l ± — -7=(vT/ 01 ± 1^02)- The energy 
splitting between two zero-energy states in vortices in a 
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spinless p x + ip y superconductor has been calculated in 
the limit of large separatio n 25 ' 26 : 



-E- 



2^ cos(kpR 



tt/4) 



(29) 



with Aq being the amplitude of the bulk superconducting 
gap, kp the Fermi momentum and £ the coherence length. 
The exponential decay of splitting is universal for all non- 
Abelian topological phase and is in fact the manifestation 
of the topological protection. 

The Berry matrix M can be evaluated and it takes the 



J 



following form: 



M 



a 

T? 



(30) 



where a = a(R) is expressed as the overlap integral of 
the two bound state wavefunctions. The form of a is not 



important apart from the fact that |a| 



-R/t 



Since both a and E± are functions of R, they are time- 
independent. We now have to solve essentially textbook 
problem of the Schrodinger equation of a spin 1/2 in a 
magnetic field, the solution of which is well-known: 



c + (T) 
c_(T) 



'cos ST- 

jug* 



i§±sin£T 



^sin£T \ 
sin £T cos S T + ift s j n £T ) 



fc+(P)\ , 



El +uj 2 \a\ 2 . 



(31) 



r 



We can translate the results into transformation of Ma- 
jorana operators: 



71 - > 

72- ^ 



^4 



£ 

-iculma 



( iwRea „ \ E + — iujIma 
^ cos £T+ - sin ST) 72 - + sin STji 

sin £T7 2 - ^ cos ST — sin £Tj 71 

(32) 



Because E + ~ A e" fl/ ^, |a| - e~ R ^ and w < A , by 
order of magnitude we can safely assume \E + \ ^> w\a\. 
In the limiting case uj — > we find 



71 — > cos £T72 — sin STji 
72 — > — sin £T72 — cos £Tji ' 



(33) 



which can be compactly written as ji — > U^iW where 

U = exp ^ - ^pj 7271 . (34) 

When we also take into account the terms containing 
w, the transformation matrix is no longer in SO(2). This 
simply means that we can no longer make the assumption 
that the system stays in the ground state subspace. 



gap. In addition, there may be low-lying bound states 
within the bulk gap, such as the CdGM states in vor- 
tices. We call them subgap states. Extended states and 
subgap bound states apparently play different roles in the 
braiding of Majorana fermions. To single out their effects 
on the braiding we consider them separately and in this 
subsection we consider excited bound states first. Since 
the energy scale involved here is the superconducting gap, 
we neglect the exponentially small energy splitting whose 
effect has been considered in the previous subsection. 

The BdG wavefunctions of excited bound states in each 
defect are denoted by — 1,2 where i labels the de- 

fects, with energy eigenvalues e\. In Appendix[B]we show 
the analytical forms of these wavefucntions for bound 
states in vortices known as CdGM bound states. If no 
other inhomogeneities are present, the wavefunctions are 
all functions of r — R^. Assuming |Rj — Rj| S> £, we 
have approximately (^xil^A'j) = % up to exponentially 
small corrections. Therefore, the BdG equation decou- 
ples for the two defects since the tunneling amplitudes 
between them are all negligible. So it is sufficient to con- 
sider one defect and we will omit the defect label i in the 
following. We write the solution to time-dependent BdG 
equation as 



B. Effects of Excited Bound States 

The concept of quantum statistics is built upon the adi- 
abatic theorem claiming that in adiabatic limit, quantum 
states evolve within the degenerate energy subspace. Go- 
ing beyond adiabatic approximation, we need to consider 
processes that can cause transitions to states outside the 
subspace which violates the very fundamental assump- 
tion of adiabatic theorem. In the case of the braiding of 
Majorana fermions in superconductors, there are always 
extended excited states in the spectrum which are sepa- 
rated from ground states by roughly the superconducting 



where we have defined as the particle-hole conjugate 
state of with energy eigenvalue —£\. We will focus 
on the minimal case where only one extra excited state 
is taken into account. Actually, in the case of bound 
states in vortices, due to the conservation of angular mo- 
mentum, zero-energy state is only coupled to one excited 
bound state (See Appendix |B|) and to the leading order 
we can neglect the couplings of the zero-energy states 
to other excited states as well as those between excited 
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states. The time-dependent BdG equation reduces to 

(36) 



ic = -/3c x + /3*c A 
ic\ = {e\ - a)c\ - P*c , 
ic\ =-{e\ + a)c x - /3c 



where we have defined the components of the Berry ma- 
trix as fix = i(^>o\i>x)-,ax = i(* a |*a)- Expressions of /3 
and a can be found in Appendix [EJ In the following we 
suppress the subscript A and shift the energy of excited 
level to eliminate ax- £\ — > ex + ax- Solving the BdG 
equation (|36|) . we find 



7^ 



4^_ :J ^. 2\/2^esin 2 ^~ \/20sw. Et 



1 — mii "t' 

2\/2e/3sin 2 ^ / 
E2 -7+ cosM 



E 2 

4/? 2 sin 2 f \ 
£ 2 I 4 



e sin 



ground state subspace and excited states as long as they 
are both localized and can be considered as a composite 
qubit. 

C. Effects of Excited Extended States 



We next consider the effect of excited states that are 
extended in space. Usually such states form a continuum. 
The scenario considered here may not be very relevant to 
Majorana fermions in 2D p x + ip y superconductor since 
bound states in vortices often dominate at low energy. 
However, for Majorana fermions in one-dimensional sys- 
tems, zero-energy state is the only subgap state and cou- 
pling between zero-energy state and continuum of excited 
states may become important. 

Again we write time-dependent BdG equation in the 
instantaneous eigenbasis 



\/2B$mEt^ esinEt? 
V ~> F 7 B — £ + cos V 



E 



E 



E 



(37) 



which should be followed up by the basis transformation 
B. Here E = y/e 2 +4(3 2 . 

By using (|24[) we can work out explicitly how the 
ground state wavefunctions transform. Physically, the 
non-adiabatic process causes transitions of quasiparticlcs 
residing on the zero-energy level to the excited levels. 
Superficially these transitions to excited states signifi- 
cantly affect the non-Abelian statistics, since the parity 
of fermion occupation in the ground state subspaces is 
changed as well as the quantum entanglement between 
various ground state o 33 ' 38 . This can also be directly seen 
from (j3"7| since starting from |g) the final state is a su- 
perposition of |g) and (ijd^lg). So we might suspect that 
errors are introduced to the gate operations. However, 
noticing that the excited states are still localized, the par- 
ity of the total fermion occupation in the ground state 
subspace and local excited states are well conserved. This 
observation allows for a redefinition of the Majorana op- 
erators to properly account for the fermion occupation in 
local excited states, as done in [39| . We therefore have 
generalized conservation of "local" fermion parity 



V12 = Hil2 II (1 2SJi 



(38) 



shared by defects 1 and 2. It is not difficult to see that 
these generalized fermion parity operators(thus their ex- 
pectation values) transform exactly according to Ivanov's 
rule, thus naturally extending non-Abelian statistics to 
an enlarged Hilbert space. From the perspective of mea- 
surement, to probe the status of a topological qubit it is 
necessary and sufficient to measure the fermion parity as 
defined in (1381) . It is practically impossible(and unneces- 
sary) to distinguish between the fermion occupations in 



(39) 



ico.i = ~y] (Pi\cx - P{\Cx) 

A 

«C ,2 = - ^ (^ACA - /?2ACa) 
A 

icx = {ex ~ ax)cx - Plx c o,i ~ /^2A c o,2 
icx = -{ex - a x )c x + /3iaCq,i + AuCo,2 



Here we still use A to label the excited states. As in 
the case of bound states, we ignore the coupling between 
excited states since these only contribute higher order 
terms to the dynamics of zero-energy states. In another 
word, we treat each excited state individually and in the 
end their contributions are summed up. 

To proceed we need to determine matrix elements ^ia 
and /?2A- We will consider the spinless one-dimensional 
p-wave superconductor as an example. The zero-energy 
states localize at two ends of the ID system which lie in 
the interval [0,L]. We also assume reflection symmetry 
with respect to x = L/2. Without worrying about the 
tunneling splitting, we can consider the two ends near 
x = and x = L independently. Then we make use of the 
fact that BdG equations near x — and x = L are related 
by a combined coordinate and gauge transformation x — > 
L — x, A{x) — > —A(L — x). Therefore, the bound state 
and the local part of the extended states near x = and 
x = L are related by gauge transformations. Based on 
these considerations, we should have /3ia = fax = Px, 
up to exponentially small corrections. Similar argument 
also applies to vortices in 2D p x + ip y superconductors. 

As mentioned above, we will make the approximation 
that each excited state can be treated independently. So 
we consider the effect of one of the excited states first 
and omit the label A temporarily. Again we write d = 
-l={i;+ifi). Without loss of generality we also assume that 
is real. We find from the solution of time-dependent 
BdG equation that 
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7i 



72 



4/3 2 2 Et 



1 



Aj3 2 2 St 



2a/2/3£ sin 2 =f ? ^2/3 sin M . 



E 2 



E 



4/3 2 2 
- |1 -^ Sm ~2~ ' 71 



4/3 2 



£» ,8m T 72 



Si „ 2 V2/3e sh 



2 _Et 



v^/JsinEt 



S 2 



2^/32 sin 2 I* N / 

^ -(-71 + 72) + I cosM + 



3 2 sm 2 ^\~ \/2esinM„ 



£-2 



-n 



(40) 



2V2/3sin 2 Bt 



„ . £ sin i?i £ _ „ 

— (-71+72) H ^ — £,+cosEtr] 



Here again i? = e 2 + 4/3 2 . 

At first glance the physics here is very similar to what 
has been discussed for local bound states: non-adiabatic 
transitions cause changes of fermion parity in the ground 
state subspace. The crucial difference between local 
bound states and a continuum of extended states is that, 
in the former case, local fermion parity is still conserved 
as long as we count fermion occupation in the excited 
states while in the latter, it is impossible to keep track 
of the number of fcrmions leaking into the continuum so 
the notion of local fermion parity breaks down. These 
non-adiabatic effects may pose additional constraints on 
manufacturing of topological qubits. Let's consider to 
what extent the braiding statistics is affected. A useful 
quantity to look at here is the expectation value of the 
fermion parity operator in the ground state subspace, 
namely (Vq) = (17172)- 

Suppose at t = we start from the ground state |g) 
with even fermion parity (g|"Po|g) = 1 an d the excited 
level is unoccupied, too. After the braiding at time T 
the expectation value of Vq becomes 

8/3 2 _ : _ 2 ET / , 2e . 2 ET 



(V (T)) = l~j^sm'—\ [ l + -sm'—). (41) 

where V {T) = U^(T)VqU{T). This confirms that 
fermion parity is not conserved anymore. For |/3| <C e, 
the coupling to excited state can be understood as a small 
perturbation. (Vq) only slightly deviates from the non- 
perturbed value. In the opposite limit \/3\ ^> e, (Vq) can 
oscillate between 1 and —1 so basically fermion parity is 
no longer well-defined. 

Now we can sum up the contributions from each ex- 
cited state and (|4"Tj) is replaced by: 



1- 



A 



E^T 



El 



■ sin 



T 2£ A .2 

1 H sin 

E\ 



E\T 



(42) 



IV. DISCUSSION AND CONCLUSION 



In conclusion, we have considered the braiding of non- 
Abelian anyons as a dynamical process and calculated 



the corrections to non-Abelian evolutions due to non- 
adiabatic effects. We discuss several sources of non- 
adiabaticity: first of all, tunneling between non-Abelian 
anyons results in splitting of the degenerate ground 
states. The Abelian dynamical phase accumulated in the 
process of braiding modifies Ivanov's rule of non-Abelian 
statistics. Since the bulk of the superconductor is fully 
gapped such corrections are exponentially small. In the 
context of TQC such deviations from Ivanov's rule are 
sources of errors in single-qubit quantum gates. 

Second, we consider dynamical transitions of Majo- 
rana fermions in the zero-energy ground states to ex- 
cited states. The effects of such non-adiabatic transitions 
strongly rely on whether these excited states are bound 
states with discrete spectrum and localized at the same 
positions with the Majorana bound states, or they extend 
through the whole bulk and form a continuum. Gen- 
erally speaking, non-adiabatic transitions mix the zero- 
energy ground states with other excited states and it is 
questionable whether the quantum entanglement crucial 
to non-Abelian statistics is still preserved. In the for- 
mer case where excited states are localized, we are still 
able to define conserved fermion parity stored in these 
low-energy bound states. Non-Abelian statistics can be 
generalized once we enlarge the Hilbert space to include 
all local bound states. In the latter case, the situation 
is dramatically different because the notion of fermion 
parity in the low-energy Hilbert space no longer makes 
sense once extended states above the bulk gap are in- 
volved. We characterize the loss of fermion parity in such 
non-adiabatic transitions by the expectation value of the 
"local fermion parity" operator. This can be viewed as 
the dissipation of topological qubit resulting from cou- 
plings to a continuum of fermionic states. We have thus 
quantified the expectation that a zero-energy Majorana 
mode will decay if it is put in contact with a continuum 
of fermionic states(e.<7., electrons). 

Although the underlying technological motivation for 
topological quantum computation is that quantum error 
correction against continuous decoherence is unnecessary 
as a matter of principle in topological systems since de- 
coherence due to local coupling to the environment is 
eliminated, other errors, such as non-adiabaticity consid- 
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ered in this work, would invariably occur in all quantum 
systems in the presence of time-dependent quantum gate 
operations. In addition, braiding is the cornerstone of the 
strange quantum statistical properties which distinguish 
non-Abelian anyons from ordinary fermions and bosons. 
Our work, involving the non-adiabatic corrections to any- 
onic braiding, is therefore relevant to all current consid- 
erations in the subject of Ising anyons whether it is in the 
context of the observation of the non-Abelian statistics 
or the implementation of topological quantum compu- 
tation. In particular, non-adiabaticity in the Majorana 
braiding in the specific context of non-Abelian topolog- 
ical superconductors as discussed in this paper, may be 
relevant to various recently proposed Majorana interfer- 
ometry experiments involving vortices in 2D^2r— . 

We now speculate about the physical sources of non- 
adiabaticity of braidings. Throughout our work we have 
focused on the intrinsic non-adiabaticity originating from 
the fact that braidings are done during a finite interval of 
time. For such effects to be appreciable, the time-scale 
of braidings has to be comparable to At ~ where 
AE is the energy gap protecting the anyonic Majorana 
modes. For the corrections from the tunneling of Majo- 
rana fermions, At is exponentially large so it is usually 
legitimate to neglect the tunneling effect. On the other 
hand, if we arc interested in corrections from states above 
the gap, then the relevant time scale is K/A. Using (l4"2")l 
one can estimate the non-adiabatic error rate in perform- 
ing gate operations. 

These considerations also apply to other types of non- 
adiabatic perturbations as long as they have non-zero 
matrix elements between zero-energy states and excited 
states. In particularly disorder scattering, which is un- 
avoidable in solid state systems, can be a source of non- 
adiabaticity 3 ®. As the non-Abelian anyons are mov- 
ing, the disorder potential seen by the anyons changes 
randomly with time so it can be modeled as a time- 
dependent noise term in the Hamiltonian which may 
cause dephasing. Other possible perturbations in solid 
state systems include collective excitations, such as 
phonons and plasmons(or phase fluctuations). We leave 
the investigation of these effect for future work. The for- 
malism developed in the current work can, in principle, 
be used to study these dephasing errors. 
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Appendix A: Proof of Eq. (gU) 

This result already appears in Ref. [{|. To be self- 
contained we give a proof here. Assume that D is not 
singular. Then there exists a real orthogonal matrix O 



such that 



ODO T 



( o e l 

-0i o 



V 



9 N 
-0 n / 



(Al) 



where 6 n > 0. We can find a new set of Majorana oper- 
ators 



b N 



( £1 \ 

C2 



O 



C2JV-1 
\ C 2 N I 



(A2) 



that allow us to write c T Be in a canonical form: 
22 c k D k id = 2 2J m b m b' m . 



(A3) 



kl 

Therefore we have 
1 



U = exp ( —c Dc 



n 



1 if 



(A4) 

The transformation of b ni b' n under U can be easily 
worked out: 



Ub n U* = cos 9 n b n - sin6>„&' 



tfb'Jj^ = sm6 n b n 
We write it in matrix form: 



cos 9 n b' n 



(A5) 



/cos 6\ — sin 6i 
sin 6*i cos 9\ 



Ubfr = 



Notice that 



cos ( 
sin f 



— sin 9 r 
cos 6„ 



exp 



cosftjv — smc/jv 
sinfljv cos^jv 







(A6) 



(A7) 



so the big matrix is essentially e 9 . Now the rest of the 
proof is rather straightforward: 



Ucfr = T e~ B b = T e e Oc 



e~ D c. 



(A8) 



Thus we have e~ D = V. 

Formula (|24[) gives the explicit form of unitary trans- 
formation corresponding to a cyclic time evolution. How- 
ever, analytically finding matrix D as the logarithm of 
the matrix V is still quite difficult. In fact, for the pur- 
pose of calculating the operator expectation values, it is 
not necessary to know explicitly D. To illustrate this 
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point, consider the ground state expectation value of 
A = Cjj Ci 2 ■ ■ ■ c- H before and after braiding. We denote 
the starting ground state by \ip). Then the expectation 
value on the transformed ground state after braiding is 
given by 

(il>\U*Au\il>) = (^c n U ■U^c l2 U---U^c H U\^), (A9) 
where WcU can be easily evaluated: 

tfcU = e D c = V T c. (A10) 

Appendix B: Bound states in vortices 

In this appendix we present the analytical solutions of 
CdGM bound states 4 ^ in spinless p x + ip y superconduc- 
tor o 4 ' 38 i 46 in the weaking-coupling limit (A <§C Ep). We 
model the p x + ip y pairing by the following gap operator: 

A(r) = --3-{A(r),a c + ^}. (Bl) 

Kp 

A vorticity 1 vortex sitting at the origin is represented 
by the following order parameter 

A(r) = A(r)e lS , (B2) 

where (r, 6) is the polar coordinate of 2D plane. 

The (unnormalized) wave function of bound states are 
given by 

Here fi(k±r) is defined to be 

f t (k ± r) = Ji(k±r)e-^K dr ' A(r '\ (B4) 



where k± — fc M ± meifk^ and the corresponding energy 
levels 

A 2 

ei = -(I-l)wo,a;o~-£. (B5) 

£F 

In weak-coupling limit we neglect the difference be- 
tween /i and Fermi energy Ep so k± rs kp. 

Zero energy state corresponds to Z = 1. Written out 
explicitly, it has the form 

*o(r) = (_£-*) Ji{kpr)e-^ fo dr ' A ( r '\ (B6 ) 

Notice that we have multiplied the wavefunction by i to 
satisfy the Majorana condition: Vi = \ = u*_ v The first 
excited state has 1 = 0: 

(u l=0 (r)\ _ f Jo(kpr) \ / -dr'A(r') , R7 ^ 

Having all these analytical expressions, we can calcu- 
late the Berry connections. We skip the algebra and the 
result is given by 

i{^i=i,i\^i,i) = -i^A[(v x -iv y )8i-(v x +iv y )5i- 2 ], (B8) 
where v = R and constant A is the following integral: 

/•oo 

A= rdr fi(k F r)d r fo(k F r)- 
Jo 

j r&rh(k F r)(d r + ^jf 2 (k F r). (B9) 
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